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Abstract 

The calculation of the polarization in ferroelectric thin films is performed 
using an analytical solution of the Euler-Lagrange differential equation with 
boundary conditions with different extrapolation lengths of positive sign on 
the surfaces. The depolarization field effect is taken into account in the model 
for a short-circuited single domain film, that is a perfect insulator. It is shown 
that the calculation of the polarization and other properties profiles and av- 
erage values can be reduced to the minimization of the free energy expressed 
as a power series of the average polarization with a renormalized coefficient 
which depends on temperature, film thickness, extrapolation lengths, and a 
coefficient for the polarization gradient term in the free energy functional, the 
depolarization field being also included into the renormalized coefficient. The 
function defining the space distribution properties is calculated as well and its 
amplitude is shown to coincide with the average polarization. The detailed 
calculations of the spontaneous polarization, dielectric susceptibility and py- 
rocoefficient is performed. The divergence of the dielectric susceptibility and 
pyrocoefficient for critical parameters of the thickness induced ferroelectric 



phase transition, namely at temperature T c i and critical length l c , is shown 
to exist with and without the depolarization field contribution, although the 
values of T c i and l c are different in both cases. The detailed analysis of the 
depolarization field space distribution and of this field dependence on tem- 
perature and film thickness is performed. 

I. INTRODUCTION 

The non homogeneity of the polarization and other properties is known to be a charac- 
teristic feature of ferroelectric thin films. The physical reason of this phenomenon is related 
to the influence of the films surface, where, e.g., the polarization differs essentially from 
that in the bulk. Because of this polarization heterogeneity, the phenomenological approach 
to its calculation was based on the solution of Euler-Lagrange and Lame type differential 
equations respectively for polarization and dielectric susceptibility with specific boundary 
conditions [1,2]. Because the analytical solution of these equations appeared to be cumber- 
some, their solution was performed numerically in most of the papers, mainly for BaTiC>3 
and PbTi0 3 (see e.g. [3,4,5]). Therefore the analytical calculations of the film properties ap- 
peared to be more complex than in bulk materials, where the minimization of the free energy 
as a power series of homogeneous polarization makes it possible to calculate the properties 
analytically by a simple way. The real situation in films is even more complex because of 
the depolarization field contribution, which, contrary to the bulk, is non zero even in short 
circuited films because of the polarization non homogeneity. Although the depolarization 
field is able to completely destroy ferroelectric polarization [6] there are only a few works 
devoted to calculations of its contribution to the film properties [7]. The majority of the au- 
thors performed the calculations without taking into account the depolarization field effect 
[3,4,5]. It might be supposed that the flattening of the polarization via the depolarization 
field contribution, as shown in [7], could create conditions for a simplification of the thin film 
properties description, particularly on the basis of a free energy expansion similar to that in 
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bulk materials, but with renormalized coefficients. Such approach was successfully applied 
recently to the description of thin films [8] and nanomaterials [9]. However, up to now, 
there was no fundamental background for such a simplified approach. In addition, nothing 
was discussed about the properties profiles, which cannot be obtained from conventional 
expansion of the free energy as a power series of the polarization. 

In the present work we show for the first time that polarization and other properties 
averaged over the film thickness can be found by the minimization of the free energy with the 
same form as known for the bulk materials, but with the coefficient before square polarization 
depending upon temperature and film thickness, and featuring the extrapolation length 
and coefficient before polarization gradient. The space distribution of the polarization is 
calculated using the analytical solution of the Euler-Lagrange equation. The amplitude of 
this distribution is shown to be equal to the average polarization with a good accuracy 
and can thus be obtained from the aforementioned free energy. The depolarization field 
contribution is taken into account in the model which was proposed in [7] for single domain 
ferroelectric films under short circuit conditions, ferroelectric treated as a perfect insulator. 
This model seems to be reasonable because the thinner the film the better the conditions 
for single domain state appearance (see e.g. [2], [10]). In addition, for many perovskite 
ferroelectrics, the conductivity can be small enough [11]. To analyze the depolarization field 
effect we perform a detailed comparison of the properties calculated with and without taking 
into account the depolarization field. 

II. BASIC EQUATIONS 

Let us consider a thin ferroelectric film polarized along the z-axis (i.e. P z ^ 0, P x = P y = 
0), which is perpendicular to the surface of the film. This type of polarization can appear 
as a result of self-polarization of a film grown under special technological condition without 
application of any external electric field [12,13]. Since the type of substrate and electrode 
was shown to be important, the mechanical strain related to the mismatch of the substrate 
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and the film lattice constants and the thermal expansion coefficients can contribute to a 
self-polarization phenomenon. 

The polarization equilibrium value can be obtained in the framework of the phenomeno- 
logical thermodynamic theory from the minimum of the functional of the free energy [14] . In 
the considered case of polarization perpendicular to the surface of the film, it is necessary to 
take into account the depolarization field, which is proportional to the value of polarization 
and has the opposite direction, in such a way that it lowers or even cancels P s . In accepted 
models, the ferroelectric is regarded as a perfect insulator under short-circuit conditions. 
Taking into account the symmetry of the considered structure (i.e. the film polarization 
depends on the coordinate z only) one can write the free energy density functional for phase 
transitions of the 1st and 2nd order (7 7^ and 7 = respectively) as follows: 



^P 2 M + faz) + lPl{z) 



S fdPJz) 



+ 2i 



Ai 



6 



+ 2tt 



dz 



-E z P z {z) 



dz + 



(1) 



j J P*(z)dz -i-J P*(z)dz 



Here / is the thickness of the film, the coefficient a depends on temperature T as a = 
a (T — T c ), T c is the temperature of the transition from paraelectric to ferroelectric phase in 
the bulk ferroelectric, E z is the external electric field, Ai j2 are extrapolation lengths. The last 
two and the previous two terms represent the depolarization field energy and surface energy 
respectively. In what follows we will consider the case of positive extrapolation lengths, i.e., 
when polarization on the surface is smaller then in the center of the film. Only in this case 
surface effects can lead to a size-driven phase transition [1]. 

The minimization of the functional (1) leads to the following equation determining the 
space distribution of polarization over the thickness of the film: 
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The last two terms in Eq. (2a) represent the depolarization field: 

Ej = -4tt \Pz-j I P z dz I . (3) 



Here the first term is the depolarization field for the free standing film, and the second 
one is the compensating field of free charges on the short-circuiting electrodes. The depolar- 
ization field (3) is of course zero when the polarization space distribution is homogeneous. 

III. CRITICAL TEMPERATURE AND THICKNESS OF THE SIZE-DRIVEN 

PHASE TRANSITION 

Equation (2a) for the polarization space distribution is a non-linear inhomogeneous in- 
tegral differential equation and is cumbersome to solve analytically. But in the paraelectric 
phase, where the spontaneous polarization P s equals zero and the polarization is only pro- 
portional to the external field, the non-linear terms in Eq. (2a) can be negligibly small. 

Therefore the solution of this linear equation can be expressed in terms of elementary 
functions. With boundary conditions (2b), the space distribution of the polarization in this 
case has the following form: 

pPE(n _ 1-^(0 p . 

= (1 + d 2 - (1 - rfi) exp(-fe)) expj-Q + (1 + di - (1 - d 2 ) exp(-fe)) exp(e - h) _ 
V?lU (l + di)(l + d 2 )-(l-di)(l-d 2 )exp(-2/i) ' U 

1 - exp(-fe) 2 + di + d 2 - (2 - d 1 - d 2 ) exp(-h) 
~ {lf} ~ h (1 + rfi)(l + d 2 ) - (1 - d 1 )(l - d 2 ) exp(-2/i) • 

Here all 'z' indexes are omitted, angle bracket means z coordinate averaging and the 
following definitions are introduced: 

The average polarization of the film can be easily derived from (4): 

" M ~ («•) 
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Note that in the case d\ = d 2 the space distribution of the polarization in the paraelectric 
phase (4) and its average value (6) coincides with those obtained earlier [7]. 

The derivative x — (dP/dE)\ E=0 represents the linear dielectric susceptibility of the film. 
The analysis of expression (6) has shown that the value of x is positive when either / < 
(i.e. a > 0, T > T c - paraelectric phase of bulk ferroelectrics) or / > and $ > / (T < T c - 
ferroelectric phase of bulk ferroelectrics). The dielectric susceptibility diverges for $ = /, 
and the considered system undergoes a phase transition. 

Because of the / and <3> dependences on thickness and temperature, equation $ = / 
determines accordingly the critical temperature T c i at fixed thickness or critical thickness l c 
at fixed temperature. 

Note that these critical parameters correspond to the points where the ferroelectric phase 
loses its stability (1st order phase transition) or to the critical temperature and thickness 
for a 2nd order ferroelectric phase transition. The parameters / and Id can be estimated 
respectively as \ je and r c j y/e where e and r c are the dielectric permittivity and correlation 
radius of the bulk ferroelectric in the paraelectric phase respectively. Therefore for conven- 
tional values of these quantities / « 1 and I /Id ~ h 3> 1. These inequalities essentially 
simplify the form of </?(£) and $ so the equation determining the phase transition point can 
be rewritten as follows: 

1 
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The critical thickness is readily obtained from this expression: 
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(9) 



I VAi + Id Ai + Id 

Here / (0) is the correlation length at zero temperature. As it follows from Eq. (9), if 
the thickness value is less than l c (0) = Zq(0)/((Ai + Id)^ 1 + (A2 + Id)^ 1 ), then the critical 
temperature becomes negative, i.e. phase transition vanishes. 
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Therefore the ferroelectric phase transition in a film can be achieved by changing the 
film thickness at some fixed temperature or by varying the temperature of the film with 
given thickness. As a matter of fact the curve described by equation (9) determines the 
phase boundaries between the paraelectric and ferroelectric phases, i.e. the phase diagram 
of the film in temperature - thickness coordinates, which is depicted in Fig.l by solid curves 
for different values of the extrapolation lengths. The critical thickness values / c (0) are the 
intersections of the curves with the abscissa axis. Namely, at / < Z c (0), the paraelectric 
phase (polarization P s = 0) exists in the whole temperature region, while at I > Z c (0) both 
ferroelectric phase (P s ^ 0) at T < T c i and paraelectric phase at T > T c i are present. It 
should be noticed that paraelectric and ferroelectric phases should coexist in the FE region 
of Fig. 1 for first order phase transitions. To clear up the effect of the depolarization field 
we also conducted calculations without any depolarization field contribution. The results of 
the calculations of T c i and l c for 2nd order phase transitions are given in the Appendix (Eqs. 
(A5), (A6)) and in Fig. 1 (dashed curves). 

It is seen from Fig. 1 that the depolarization field decreases the critical temperature and 
increases the critical thickness. These effects are stronger when decreasing the extrapolation 
lengths values (compare solid and dashed curves in Fig.l). It is clearly seen, from the 
inset of Fig. 1, that (T c — T c i) ~ 1/7 (solid curves), while without the depolarization field 
contribution (T c — T d ) ~ l/l 2 (dashed curves). A more detailed critical thickness dependence 
on the extrapolation lengths is represented in Fig. 2. One can see that the greater the di 
and g?2 values, the smaller the l c value. The latter can be related to the flattening of the 
polarization space distribution on increasing the values of the extrapolation lengths, that 
results into a critical thickness decrease and in a critical temperature increase (compare 
curves 1, 2 and 3 in Fig.l). 

When inequalities $ < / and / > are valid, the polarization derived by equation (4) 
becomes negative, and the non-linearity in equation (2a) cannot be neglected. 
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IV. THE FREE ENERGY 



P™(0 = P(1 -¥>(£)), <p(0 



In the ferroelectric phase (the region I > Z c (0) and T <T d ) the non- linearity should be 
taken into consideration. The simplest way to take it into account is the direct variational 
method. We choose solution (4) as a trial function and an amplitude factor will be treated 
as a variational parameter. The condition of FE existence h > (Ui + U 2 )/ f can be fulfilled 
because of small / and l d values. Allowing for h 3> 1 we will look for a polarization space 
distribution in the ferroelectric phase with the following form: 

C/iexp(-0, 0<£< h/2; 

0, Z~h/2- (10) 

C/ 2 exp(£-/i), 0</i-f«/i/ 2 - 
Here P is the variational parameter that represents the amplitude of the polarization 
space distribution. 

After calculation of integral (1) with the trial function (10), one can easily obtain the 
free energy density as follows: 

F = a (l - ^ (1 - A,)^ + (3(1 - B)^ + 7 (1 - C)^ - EP{1 - A,) (11a) 

Keeping in mind that the polarization average value is (P) = P(l — Ai)one can rewrite 
Eq.(lla) as: 

(l-A/V jP) 2 (l-B) (P) 4 1-C (Pf 

= a ^A, — T + ^(T^^^ + 7 JTTa^— ~ e{p) - (llb) 

where the following definitions have been used 

= U l + U 2 = G{Ui) + G{U 2 ) = Qm + Qm (u) 
1 h ' 12h 60h 1 ' 

G{Ui) = A8Ui - 36U? + 16U? - 3Uf, 

Q(Ui) = 3Q0Ui - + AOOUf - 225E/? + 72C/f - lOC/f 

It is worth to underline that the terms in Eq. (1) which correspond to surface energy, 
depolarizing field and polarization gradient contribute to the first term in Eqs. (11). There- 
fore the main peculiarities of thin film properties have to be related to the coefficient of the 
second power of the polarization. 



It can be seen that Eqs. (11a) and (lib) have the form of power series of the amplitude 
of the polarization space distribution and average polarization respectively. These equations 
can be easily rewritten in the conventional for bulk ferroelectric form 

p2 p4 p6 

F = a h b h c E'P (13a) 

2 4 6 

and 

F = aj&- + + Cl { -^-E(P), (13b) 

where the coefficients a, b, c and a±, b±, c\ can be obtained from the comparison of Eqs. 
(lla,b) and (13a,b), namely 

a = a(l-^J(l-A 1 ), b = P(l-B), c = 7 (l-C), E' = E(l - A,) (14a) 

and 

a(l-A 1 /f) 13(1 -B) 7 (1 -C) 

ai= ' bl = ^A^> Cl= o^F (14b) 

In the general case, the coefficients (14a,b) depend on the temperature, the film thickness, 
the extrapolation lengths and the coefficient before polarization gradient, as it follows from 
Eqs. (12), (7) and (5). At first glance these dependencies seem to be rather complex, 
although they have been already simplified due to inequalities / = 1/e <ti 1, h — l/l d ^$> 1. 
A further simplification of the coefficients (14) was made by rewriting them via the critical 
temperature T c i = T c — (U\ + L^/cko • 4^/^// or via the critical thickness l c /ld = {Ui + 
U 2 )4:ir I (a (T c — T)) (see the left hand side part of Eq. (8)). The first or the second type of 
the coefficients representation can be useful when studying respectively the temperature (at 
given thickness) or thickness (at fixed temperature) dependencies of the film properties. In 
the first case the coefficients a in (14a) and ai in (14b) can be rewritten as: 



a = a (T - T d ) 

and 



1~^{T C -T d ) 



(15a) 
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a x =a (T - T cl ) 



1 



(15b) 



1-M(T C -T d )' 



An estimation of the second multipliers in Eqs. (15a) and (15b) has shown that they are 
very close to unity (for ferroelectric phase transition of displacement type a /(47r) < 1CT 5 
with T c - T cl < 10 2 and for order-disorder type a /(4vr) « 10~ 3 with T c - T cl < 10 2 ) and 
so 1 — Ai pa 1 everywhere in Eqs. (14a) and (14b). Similar estimations have shown that 
Bel and C<lso that with a good accuracy 



When studying the thickness dependence of the film properties at fixed temperature, 
substitution of critical thickness Z c (8) into Eq. (14) yields: 



the coefficients b, bi, c, c\ and E' being the same as in (16a). We should underline once 
more that the critical temperature T d depends on I and the critical thickness l c depends on 
T, and both of them are functions of the extrapolation lengths and parameter S (see Figs. 
1, 2 and Eqs. (8), (9)), the contribution of depolarizing field also being included. Because 
of that, all the physical properties (polarization, dielectric susceptibility, pyrocoefficient, 
entropy, thermal capacity etc.), which can be obtained from conventional minimization of 
Eqs. (13a), (13b), have to be dependent on the aforementioned characteristics. 

It should be also underlined, that Eqs. (13a), (13b) and (16a), (16b) are valid both in 
paraelectric (T > T c i, I < l c ) and ferroelectric (T < T c i, I > l c ) phases. This is because we 
chose the trial function for the description of the polarization profile in the form obtained 
for the paraelectric phase (see Eq. (4)) and used the inequality /i » 1, / « 1, which are 
valid practically for any thin film, as discussed earlier. 

Since the formulas which express the physical properties via the free energy coefficients 
for both first and second order phase transitions are of the common knowledge (see e.g. [15]), 
we shall consider only polarization, dielectric susceptibility and pyrocoefficients with special 



a = ai = a (T - T d ), b = b 1 = (3, c = c x = 7, E' = E. 



(16a) 




(16b) 
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attention to their thickness and temperature dependencies, which are the characteristic 
features of ferroelectric thin films. For sake of simplicity, we shall represent, in what follows, 
these dependencies for the second order phase transition. Special attention will be payed 
to depolarization field contribution to find out the cases and conditions at which these field 
effects are more or less important qualititavely or quantitatively. 

V. THE SPECIFIC BEHAVIOR OF THE PHYSICAL PROPERTIES OF 

FERROELECTRIC THIN FILMS 

A. The space distributions of the physical properties 

The profiles of the physical properties are defined by the function (1 — </?(£)) in Eq. (10) 
and their amplitudes can be found by conventional minimization of free energy (13a) with 
the coefficients (16a), (16b). More particularly, for second order phase transitions at E — 0: 



/ ao(T-T d ) 1 a _ 

V P ' X 2a (T d -TY (32P V3V7£=T' 



P iVl il /c) ' X 2Kl d /l c (l-l c /iy U 2^Ki a Jl c {l-lJl) ^ 
for the ferroelectric phase (T < T d , I > l c ) and 

^°'^;MrW n ^ (18a) 

^°-^ A WA-i) ' n ^° (18b) 

for the paraelectric phase (T > T ch I < l c ). 

Here P, x an d n are the amplitudes of the polarization, dielectric susceptibility and 
pyrocoefficient which have to be multiplied by (1 — </?(£)) to obtain these properties pro- 
files. The formulas (17a), (18a) and (17b), (18b) can be applied when studying respectively 
temperature (at fixed film thickness) and thickness (at fixed temperature) dependencies. 
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The space distributions of these quantities are depicted in Figs. 3, 4, 5 for several film 
thickness at fixed low temperature, P s0 , Xo, n being the quantities value at I — > oo (see 
Eq. (17b)). One can see that the amplitude of the susceptibility increases much faster 
than that of the pyrocoefncient at I — > l c while the amplitude of spontaneous polarization 
decreases with the decrease of the film thickness, the shape of all these quantities profiles 
being completely the same. It should be underlined that the profiles are asymmetrical for 
different extrapolation lengths on the two surfaces. It can be easily seen that the profiles 
remain smooth in the most part of the film even for thin films with thickness close to the 
critical value l c (see solid curves 1 and 2). This phenomenon can be explained by the 
influence of the depolarization field that tends to flatten the polarization space distribution 
as it will be shown later. 



The thickness and temperature dependencies of the properties can be obtained directly 
by minimization of the free energy (13b) with respect to Eqs. (16a), (16b) and accordingly, 
they can be calculated on the basis of Eqs. (17) and (18) for the ferroelectric and para- 
electric phases respectively because the average values coincide with the properties profiles 
amplitude. As an illustration, we depicted in Figs. 6 and 7 the temperature dependence of 
the inverse dielectric susceptibility and pyrocoefncient, the latter being proportional to the 



spontaneous polarization (see Eq. (17a)). In Figs. 6, 7 we introduced P s q(0) = Ja T c /f3, 



Xo(0) = l/(a T c ) and n = y/dio/ (2^ f3T c ) which are these properties at T = in the bulk 
material. The thickness dependencies are depicted in the insets to Figs. 3, 4, 5. In particu- 
lar, the slope of the strait line in the inset to Fig. 4 gives the value of the parameter 2Kl d /l c 
that defined both susceptibility, pyrocoeflicients and polarization (see Eq. (17b) and inset 
to Fig. 5). 

A characteristic feature of the properties temperature dependencies is the shift of the 
average properties anomalies to lower temperature on decreasing the film thickness. Since 



B. The average values of the physical properties 
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the divergences of x an d n are related to the temperature of the thickness induced ferro- 
electric phase transition (see Eqs. (17), (18)), the aforementioned shift is related to the T d 
dependence on the film thickness (see Fig. 1). The shift of T d due to the depolarization 
field effect (compare dashed and solid lines in Fig. 1) defines this field effect on the average 
properties of the ferroelectric films. Since the intersections of the curves with abscissa axis 
in Fig. 1 defines the / c (0)/7o(0) value, the shifts of / c (0)//o(0) to larger values due to the 
depolarization field contribution is clearly seen. 

The estimations of Z (0) = ^j5/(a T c ) and l d = ^5 /(Ait) values for the parameters of 
PbTi0 3 and BaTi0 3 [3,5] gave respectively l (0) = 15 A, l d = 0, 6 A and / (0) = 40 A, 

o 

l d = 2 A which makes it possible to obtain the / c (0) value from the intersections in Fig. 1 
or by the calculation of l c on the basis of Eq. (8). For instance, for extrapolation lengths 

o o 

Ai = A 2 = /o(0) / c (0) = 30 A and 80 A with depolarization field contribution and without 

o o 

it Z c (0) = 24 A and 63 A for PbTi0 3 and BaTi0 3 respectively. Critical thickness has to 
increase with decrease of Ai and A 2 values and at T ^ namely l c (T) = / c (0)/(l — T/T c ) (see 
Eq. (8)). Because of the scattered 5 values given in different works, the obtained l c values 
should be considered as estimations only. But it is obvious that the larger T c the smaller l c . 
An estimation of the T d value from Fig. 1, e.g. for ///o(0) = 10 and parameters of curve 

o 

2, gives T c \ = 0, 8T C for the films with thickness 200-500 A. Therefore for such thin films 
the shift AT C = T c - T cl = 0, 2T C and, e.g., for BaTi0 3 and PbTi0 3 , AT C = 24° C and 
AT C = 108° C respectively, which means that T d can be rather close to the T c value. For 
thicker films T d is expected to be even closer to T c (see Fig. 1). 

C. The depolarization field and its influence on the physical properties 

The depolarization field is given by Eq. (3), which after substitutions P z = P(l — </?(£)) 
and (P) — P(l — Ai) (see section 4) can be rewritten as 

E d M) = -^P{A 1 -^)). (19) 
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Although A 1 value is small, it cannot be neglected in the region where </?(£) = (see Eq. 
(10)). In particular in the very center of the film: 

E d = E d z {£ = h/2) = —AnPAi. (20a) 

Allowing for A\ = ld{U\ + U 2 )/l one can rewrite Eq. (20a) with respect to (17b) as: 



E d = -K ld 



Kl *( l '--?\ (20b) 



l\ (3l c \l 



or 



K = - T)a (T c - T cl ). (20c) 

Equations (20b) and (20c) should be applied when considering E d thickness or temper- 
ature dependencies respectively. 

It is clearly seen that E d = at T = T d for arbitrary I and I = l c and I — > oo for arbitrary 
T. Between I = l c and I — > oo, the absolute value of E d achieves a maximum at l m = |/ c . 



The value of depolarization field at this point is equal to = j^ a y~ a /(3 — E c0 , where 
E c q is the coercive field of the bulk material. Near the surfaces, when </?(£) ^ (see Eq. 
(10)), the small quantity A 1 can be neglected and E d z becomes thus positive, its profile being 
defined by the </?(£) dependence. The point £ cr at which the change of depolarization field 
sign takes place can be obtained easily from Eq. (19), with E d (^ cr ) = at y?(£cr) = A 1 . With 
respect to Eq. (10) this gives {l J rU 2 /U 1 )l d /l = exp(-£ cr ) or -£ cr = \n(l d /l) +]n(l + U 2 /U 1 ). 
Keeping in mind that l/ld — h >> 1, and that U\ and U 2 are quantities of the same order 
of magnitude (e.g. U 2 jU\ = 1 at Ai = A2) one can conclude that for symmetrical boundary 
conditions 

&r = ln^. (21) 

Therefore at C, < £ cr E d (^) > 0, while < at £ > £ cr . Similar conditions on the right 

hand side of the film can be obtained by substitution of (h — £ cr ) for £ cr in Eq. (21). 

To illustrate the E d behaviour we have represented its profile and the thickness depen- 
dence of E d (z = 1/2) in the inset b) of Fig. 8. One can see that, because there is no 

14 



depolarization field in the film center at critical thickness, increases on increasing l/l c , 
achieves maximum at l c /l m = | and at I — > oo E d z — > again. More details for the region 
Ez(£) < are given in inset a). 

The effect of the depolarization field is clearly shown in Figs. 3 to 5. The dashed lines 
are the profiles calculated without any depolarization field contribution following the way 
described in Appendix for the thinnest (curve 1) and the thickest (curve 5) films. It can be 
seen that this field influence depends on the film thickness: the thinner the film the stronger 
the change in the profile shape. Comparison of solid and dashed curves in Figs. 3, 4, 5 
shows that the depolarization field increases x an d n values, but decreases the P value for 
the thin film close to the critical thickness, while for thick films, a decrease in x an d n 
and an increase in P s are observed in most part of the film (see curve 5). It is worth to 
notice that the dashed curves in Figs. 4, 5 display small maxima near the surfaces contrary 
to solid curves which remains flat. It is clear that the influence of the depolarization field 
on the polarization is directly related to the peculiarities of depolarisation field profile (see 
Fig. 8). In particular the positive sign of -Kf(£) in the tiny regions close to the surfaces 
and negative sign outside these regions explains the increase of the spontaneous polarization 
near the surfaces and its decrease in most part of the film, leading to the flattening of the 
polarization (see Fig. 3). The difference in the behaviors induced by the depolarization field 
effect on P s on one hand, and x an d n, which are the derivative of P s , on the other hand, 
is not surprising, since, e.g., II ~ 1/P S , and an increase in P s should result in a II decrease 
and vice versa. 



VI. DISCUSSION AND CONCLUSION 

Let us discuss briefly the accuracy of the calculations and the physical model we used 
here. The variational method used for the calculation of the polarization in the ferroelectric 
phase is an approximation, its accuracy being usually a few percent. To prove this we 
performed the calculations in two important cases by other exact ways: in the vicinity of 
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the phase transition, where P s tends to zero and can be considered as a small parameter, 
and for the thick film limit, when the polarization space distribution is almost uniform. The 
solutions obtained by these methods are reported in the inset of Figs. 3, 4, 5 by dashed 
and dotted-dashed lines. It can be seen that these lines coincide with those obtained by the 
variational method with a good accuracy. 

In our calculations we did not take into account the possible shift of T c by mechanical 
strains originating from the misfit between the substrate and the film lattice constants, 
thermal expansion coefficients and growth imperfections. Because of the strain relaxation 
due to misfit dislocations, the strains are not homogeneous and in the general case their 
value are small in most part of the film [16]. We accordingly neglected the renormalization 
of T c value by mechanical strains. 

Let us discuss briefly the depolarization field effect. As it follows from the above consid- 
erations, the depolarization field does not change the general form of P s , x an d II average 
values dependences on the film thickness and temperature (see Figs. 6, 7). The depolariza- 
tion field smoothes the space distributions of these quantities, and in particular cancels the 
maxima near the film surfaces in the space distributions of x an d II obtained in the case 
where the depolarization field can be neglected (see Figs. 4, 5). 

The main effect of the depolarization field is the shift of the boundary between the 
paraelectric and ferroelectric phases on the temperature - film thickness phase diagram to 
larger thickness, as can be seen from Fig. I. This shift can be especially large for thin 
films because the critical temperature T c i obeys a T c — T c \ ~ l/l dependence when taking 
into account depolarization field contribution (see Eq. (9) and solid lines in Fig. 1) but a 
T c — T d ~ l/l 2 without this contribution (see Eqs. (A6a), (A6b) and dashed lines in Fig. 
1). Because of that a scattering of experimental critical values can be expected for films 
produced in different technological conditions and from materials with different degrees of 
purity. These factors might influence the domain structure and conductivity of the film and 
consequently the depolarization field. 

Obviously depolarization field effects in thick films under short-circuit conditions can 
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be neglected in most part of the film except close to the surfaces. But another important 
case where the contribution of the depolarization field can be neglected is related to the 
conditions on the film surfaces. More particularly, for extrapolation lengths much larger 
than the correlation length l it follows from Eq. ( A6c) that T c — T d ~ 1/7 and Eq. (A6c) 
looks like Eq. (9). 

Although we considered spontaneous polarization and other properties without Any ex- 
ternal electric field application, the case of E ^ can be considered on the basis of Eq. (13), 
i.e., in the same way as nonlinear effects in the bulk. It should be mentioned that a non 
linear dielectric response was measured in BaojSro^TiOs thin films with thickness 24-160 
nm [8]. Authors of this work explained the observed temperature and E- field dependencies 
using a conventional power series expansion of the free energy. They came to the conclusion 
that only the first term of this power series varies significantly with the film thickness or the 
temperature. This conclusion is in very good agreement with the results obtained in our 
work. The spatial distribution of the pyroelectric coefficient measured by LIMM method 
in homogeneously poled ferroelectric polymer films [17] looks like that depicted in Fig. 5. 
This may give some support to the relevance of the calculated profiles of polarization and 
dielectric susceptibility. 

The obtained results pave the way to calculations of ferroelectric film properties using a 
conventional minimization of the free energy for bulk material, but with renormalized coef- 
ficient before squared polarization. This renormalized coefficient includes the contribution 
of polarization gradient, depolarization field, polarization at the surface and the character- 
istics of the bulk material. Because this coefficient has been shown to depend on the critical 
parameters of the thickness induced ferroelectric phase transition T c i and l c in the same way 
as this coefficient temperature dependence in the bulk, namely, for the films, a = a^{T — T c i) 
or a = Kl d /l c {l c /l — 1) (see Eq. (16)), this makes it possible to obtain these critical parame- 
ters through the observation of thermodynamic properties anomalies, e.g. by measurements 
of the average dielectric susceptibility maximum position at different temperatures at fixed 
thickness or at fixed temperature and different film thickness respectively. The novel ap- 
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proach developed here is suitable for the calculation of the profiles and average values of 
the polarization, dielectric susceptibility, pyrocoefficient, entropy, specific heat etc. In or- 
der to calculate the electromechanical properties, we have to include into the free energy 
functional the interaction of the non homogeneous polarization with external strains. This 
consideration is in progress now. 



VII. APPENDIX 



In order to discuss the influence of the depolarization field, let us consider briefly the 
model without taking into account that field. In this case the equation giving the space 
distribution of polarization as a function of the thickness of the film has the form Eq. (2a) 
without the last two terms. This equation has a nontrivial solution when the electric field 
is zero. In the case of boundary conditions (2b) with positive extrapolation lengths, this 
solution exists only under the condition a < (see, e.g. |1|) and gives the following space 
distribution of the spontaneous polarization for the phase transition of the second order: 

z + z 



PAz) 



P. 



2m 



-sn 



1 + m \ l Q y/l + m 



m 



(Al) 



Here Iq= y—X/a is the correlation length, sn(u\m) is the elliptic sine function JT8J, and 
the constants m and zq have to be determined from the boundary conditions (2b). Using 
the properties of the elliptic functions it is easy to obtain equations which these constants 
satisfy: 



2(1 



/o\A + m -P 1 (arcsin(/i), m) 



(A2) 



I = £qv1 + m (2K(m) — F(arcsin(/i),m) — F(arcsin(/ 2 ), rn)) 



(A3) 



where K(m) and F((j),m) are complete and incomplete elliptic integrals of the first kind 



respectively |18| , and the following definition is introduced: 
f i = f(m,\)= — 1 + I — I - II 



\ 



2m 



V 



\ 



x, 
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{1+mY 



(A4) 
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The parameter m varies from to 1 in Eqs. (Al) - (A4), otherwise the spontaneous polar- 
ization would be a complex number. As it follows from Eq. (A3) with the condition m — > 0, 
thickness /tends to some critical value l c which corresponds to a zero value of the sponta- 
neous polarization (Al). When the thickness of the film is less than critical, the spontaneous 
polarization vanishes. The critical thickness for the considered case Ai > 0, A 2 > has the 
form: 

lc = k(^~ arctg ^j^j - arctg (A5) 

The critical temperature T c \ for a film with thickness I can be found from the condition 
/ = l c (T = T d ), which results in a transcendental equation for the dependence of T d on the 
film parameters. This equation can be simplified in several limiting cases. When the phase 
transition temperature is much smaller than in the bulk ferroelectric (T d <C T c ), i.e., for 
very thin films, the following dependence can be obtained: 



X(0) + / (0) {gi + g 2 )' 



2 N 



T « = T 4 i - {ntmf^r ) 1 ■ ' - m ^ wo) (A6a) 

Here/ c (0) and/ (0) = ^\/(a T c ) are the critical thickness and correlation length respectively 



at zero temperature T = 0, and = Aj/ \JIq(0) + Af . 

When the phase transition temperature is of the same order of magnitude as in the bulk 
ferroelectric (T d w T c ), i.e. for films with large thickness, it is easy to get from Eq.(A5) the 
following expression: 

r ^4~(^y 2 )' i>>wo) <A6b) 

Note that this expression is also valid in the whole temperature and thickness ranges for 
extrapolation lengths much smaller then the correlation length (Aj Zo)- in the opposite 
case, i.e. Aj ^> l , the critical temperature is expressed as: 

Td = Tc ( 1 "^(^ + i))' Ai>>/ ° (A6C) 
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Since a = ao(T — T c ), where T c is the temperature of the ferroelectric phase transition of 
the thick film, the nonhomogeneous polarization of the film P{z) should depend on tempera- 
ture T and external electric field E. This allows to calculate the non homogeneous pyrocoef- 
ficient n(z) = (dP (z) / dT) e=q and the linear dielectric susceptibility x{ z ) = (dP(z)/dE)E=o- 
Differentiation of Eqs. (2a) and (2b) gives the following differential equations for the pyro- 
coefficient and dielectric susceptibility calculations: 



(a + 3/3P s 2 )II - 5^ + a P s = 



dU 




n 




dU 


n 


dz 


2=0 




z=0 


dz 





z=l 



(A7a) 
(A7b) 



{a + WP 2 s )x- ^ -1 = 0; 



dx 




X 
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dz 
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Ai 


2=0 


dz 


z=l 


A 2 



(A8a) 



(A8b) 



z=l 



The general solution of Eq. (A8a), which is a particular case of the Lame equation, is the 



following [11J: 



Il(x)=CZ y y 1 (x)+C™y 2 (x)+y™(x) 
where the function yi(x) has the following form: 

yi(x) = cn(x, m)dn(x, m), 



(A9) 



(AlOa) 



Here cn{u\m) and dn{u\m) are the elliptic cosine and delta amplitude functions respectively 



y 2 (x) = [x-- E(am(x),m)) ^^+y {x) — ^ — x2 J y0K , (AlOb) 
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Here IIo is the thick film pyroelectric coefficient and the following definitions are introduced: 

z + z 

x = - /== , Vo\x) = sn(x,m). 
/ovl + m 

In Eqs.(AlOb, c) E(<p,m) and am(x) are an incomplete elliptic integral of the second kind 
and an elliptic amplitude function respectively |L8| . The constants Cf y and C% y are omitted 
because of their cumbersome expression. 

The dielectric susceptibility can be obtained in a similar way: 

X (x) = C?yi(x) + C?y 2 (x) + yf{x) (All) 

where 

rh, x /. . 1 + m — 2my (x) , »,„s 

(x) = -2 Xo (l + m) (1 _ m)2 (A12) 

Here Xo = — l/2a is the thick film susceptibility. Coefficients and differ from Cf y 
and Cf y in Eq.(14) because of the difference between y^{x) and yf l (x). 



21 



REFERENCES 

[1] Y. Ishibashi, H. Orihara and D. R. Tilley, J. Phys. Soc. Jap., 67, 3292 (1998). 

[2] Y.G.Wang, W.L.Zhong, P.L.Zhang, Phys. Rev. B, 51, 5311 (1995). 

[3] B. D. Qu, W. L. Zhong, R.H. Prince, Phys. Rev. B, 55,11218 (1997). 

[4] Y.G.Wang, W.L.Zhong, P.L.Zhang, Phys. Rev. B, 53, 11439 (1996). 

[5] G.L.Wang, S.R.Smith, J.Phys.:Cond.Matter, 7, 7163 (1995). 

[6] D.R.Tilley, Ferroelectric Thin Films (Gordon and Breach, Amsterdam, 1996). 

[7] R. Kretschmer and K. Binder, Phys. Rev. B, 20, 1065 (1979). 

[8] C.Basceri, S.K.Streiffer, A.I.Kingon, RWaser, J. Appl. Phys., 85, 2497 (1997). 

[9] B.Jiang, L.A.Bursill, Phys. Rev. B, 60, 9978 (1999). 

[10] A.M. Bratkovsky, A.P Levanyuk, Phys. Rev. Lett. 84, 3177 (2000). 

[11] G.Suchaneck, T.H.Sander, R.Kohler, G.Gerlach, Integrated Ferroelectrics, 27, 127 
(1999). 

[12] G.E.Pike, W.L.Warren, D.Dimos, B.A.Tuttle, R.Rames, J.Lee, V.G.Keramidas, 
J.T.Evans, Appl. Phys. Lett. 66, 484 (1995). 

[13] R.Bruchhaus, D.Pitzer, R.Primig, M.Schreiter, W.Wersing, Integrated Ferroelectrics, 
25, 1 (1999). 

[14] V.M. Fridkin, Ferroelectrics semiconductors, N.-Y.: Consult Bureau (1980). 

[15] M.E. Lines, A.M. Glass, Principles and Application of Ferroelectrics and related Materi- 
als, Clarendon Press, Oxford (1977). 

[16] J. S. Speck, W.Pompe, J. Appl. Phys., 76, 466 (1994). 

[17] B.Ploss, R.Emmerich, S.Bauer, J. Appl. Phys., 72, 5363 (1992). 

22 



[18] Gradshteyn I.S. and Ryzhik I.M. Tables of integrals. -Moscow: Physmatgis; 1963. 

[19] Erdelyi A. et al. Higher Transcendental Functions. -N.Y.: McGraw-Hill Book CO; 1953, 
V. 3. 



23 



FIGURES 

FIG. 1. The dependence of the critical temperature T c i on the film thickness (logarithmic scale) 
for /(0) = aoT c /4ir = 0.01 and for different extrapolation lengths values: (Xi/ld, M/ld) = (0.5, 0.6), 
(5, 3), (10, 6) for lines 1, 2 and 3 respectively with (solid lines) and without (dashed lines) the 
depolarization field contribution. The behavior as a function of the reciprocal film thickness is 
given in the inset. 

FIG. 2. The dimensionless critical thickness fl c /ld m the (di,^) space. The dj's are the 
dimensionless extrapolation lengths. The critical thickness is given as numbers close to the relevant 
plots. 

FIG. 3. The space distribution of the spontaneous polarization P s as a function of the z coordi- 
nate for the following values of parameters: / = 0.01, Xi/ld = 5, \2/ld = 3, and for different values 
of the film thickness l/l c = 1.01, 1.04, 1.16, 1.7, 4 (solid curves 1-5 respectively). The dashed 
curves 1 and 5 represent the space distribution of P s calculated without taking into account the 
depolarization field. The average value of P s 2 vs the reciprocal film thickness / is represented in 
the inset. The solid line represents the dependence obtained using the variational method, dashed 
and dash-dotted curves are obtained by expansions close to the phase transition and for the thick 
films respectively. 

FIG. 4. The space distribution of the dielectric susceptibility in the ferroelectric phase as a 
function of the z coordinate for the same values of parameters as those in Fig. 3. The dashed 
curves 1 and 5 represent the dielectric susceptibility space distribution calculated without taking 
into account the depolarization field. The inset gives the dependence of the reciprocal averaged 
dielectric susceptibility as a function of the reciprocal thickness of the film. The meaning of the 
solid, dashed and dash-dotted curves is the same as in Fig.3. 
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FIG. 5. The space distribution of II as a function of the z coordinate for the same values of 
parameters as those in Fig. 3. The dashed curves 1 and 5 represent the II space distribution 
calculated without taking into account the depolarization field. The inset gives the dependence of 
the reciprocal averaged II as a function of the reciprocal thickness of the film. The meaning of the 
solid, dashed and dash-dotted curves is the same as in Fig.3. 

FIG. 6. The temperature dependence of the reciprocal averaged dielectric susceptibility for the 
following values of parameters: /(0) = aoT c /4m = 0.01, Xi/ld = 5, M/ld = 3, and for different 
values of the film thickness l/l c (0) = 1.11, 1.43, 2, 3.33, 10 (curves 1-5 respectively). The dashed 
line depicts the temperature dependence of the bulk ferroelectric reciprocal dielectric susceptibility 
Xo(T)- The temperature dependence of the average dielectric susceptibility for the same values of 
parameters is given in the inset, where the dashed line represents the bulk ferroelectric dielectric 
susceptibility. 

FIG. 7. The temperature dependence of the averaged spontaneous polarization for the same 
values of parameters as those in Fig. 6. The temperature dependence of the bulk ferroelectric 
spontaneous polarization T/T c is plotted as a dashed line. The temperature dependence of the 
averaged II for the same values of parameters is given in the inset, where the dashed line represents 
the bulk ferroelectric pyroelectric coefficient IIo = dP s o/dT = Ho(0)^/l — T/T c . 

FIG. 8. The space distribution of the depolarization field close to the film surface for the 
following values of parameters: / = 0.01, Xi/ld = 5, \2/ld = 3, and for different values of the 
film thickness l/l c = 1.01, 1.1, 1.5, 4 (curves 1-4 respectively). The region of the plot where the 
depolarization field is negative is depicted in the inset (a). The depolarization field in the center of 
the film as a function of the reciprocal film thickness is depicted in the inset (b) for the parameters 
values given above. 
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